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Abstract 

We study ground states of two- component condensates in a harmonic trap. We prove 
that in the strongly coupled and weakly interacting regime, the two components segregate 
while a symmetry breaking occurs. More precisely, we show that when the intercomponent 
coupling strength is very large and both intracomponent coupling strengths are small, 
each component is close to the positive or the negative part of a second eigenfunction 
of the harmonic oscillator in R^. As a result, the supports of the components approach 
complementary half-spaces, and they are not radially symmetric. 

1 Introduction 

A two-component Bose-Einstein condensate (BEG) is described in terms of two wave functions 
u and V, respectively representing the first and the second component. The energy of a 
trapped two-dimensional two-component BEC is given for g = {gi, g2, du) by 

£g{u,v) = \j^^ {|Vn|2 + + + {|V^;|2 + + i^abl^} + gi2\u\M'^ . 

Here 512 is the intercomponent coupling strength, and gi (respectively 52) is the intracom- 
ponent coupling strength of the first (respectively second) component. We are interested in 
the properties of the ground state of £g when gi2 goes to infinity. We assume that both 
components have repulsive internal interactions, so 

ffi>0 and (72 >0, (1.1) 
and that the trapping potential is the harmonic function 

V{x)=uj'^\x\'^ , uj>Q. (1.2) 

We consider then {ug,Vg) as a minimizer of £g over the class of finite energy pairs with 
constrained mass 

X = ^iu,v) ; u,v £ Hy(R'^;C) , ||'u||i2(]R2) = 1, ||v||l2(]r2) = l| , 

where 



The pair {ug,Vg) satisfies the system of coupled Gross-Pitaevskii equations 




—Au + Vu + gi\u\'^u + gi2\v\'^u = Xu 
—Av + Vv + g2\v\'^v + gi2\u\'^v = fj,v 



(1.3) 



in R^, where A and fi are the Lagrange multipliers associated with the mass constraints. 

The system (1.3) is a particular case of nonlinear elliptic systems with competition. This 
denomination is due to the gi2\v\'^u and gi2\u\'^v terms, which depending on the value of gi2, 
favor solutions that coexist or that spatially separate. The segregation problem consists in 
studying the properties of the solutions when gi2 is positive and very large, since in this case, 
solutions tend to have disjoint supports. 

The segregation of two-component BECs has been experimentally observed. The achievement 
of two-component BECs constitutes an important research subject in experimental physics, 
see for example [18], [26] or [28]. Two-component BECs have been realized in several ways, 
using different types of particles: two different isotopes of the same atom, isotopes of two 
different atoms, or a single isotope in two different hyperfine states. The choice of the number 
of particles and the types of isotopes determine the value of gi and g2- The value of gi2 can be 
altered by changing the hyperfine state of a portion of the particles, via applied magnetic or 
optical fields. This allows large and tunable values for gi2 ([4], [28], [31]). Therefore, manipu- 
lating the values of the coupling strengths, it is possible to realize segregated two-component 
BECs with different spatial configurations ([18], [28]). We refer to [21] for the description of 
the physics. 

In numerical analysis ([21], [22], [25]), several patterns have been observed for the supports 
of two-component BECs. In [25], Mason and Aftalion perform a numerical analysis on the 
solutions of (1.3). They classify all the possible configurations depending on gi2 and gi / g2- 
For the segregated case, they exhibit two main types of patterns: symmetry preserving and 
symmetry breaking ones. In the first case, one component is a disc centered in the minimum 
of V, and the other component is an annulus surrounding the disc. In the second case, both 
components are close to half balls. 

In this paper, we address both segregation and symmetry breaking for two-component BECs 
in all M^. Our goal is to prove that when gi2 goes to infinity, the supports of Ug and Vg are 
disjoint; and that when gi and g2 go to zero, they break the symmetry of the harmonic po- 
tential by approaching half-spaces. To our knowledge, the are no results in the mathematical 
literature about symmetry breaking for two-component BECs in the repulsive case. 

Our method consists in two steps. First, we consider sequences = (ff" , 92 ' 512) such that gi2 
goes to infinity and g^ (respectively 52 ) converges to some nonnegative limit gi (respectively 
52)- We prove that the associated sequence of minimizers {ugn,Vg-n.) converges to a limiting 
pair, which minimizes the addition of the Gross-Pitaevskii energy of each component 




'} (1.4) 



over 
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Y = ^{u,v) £X;u-v = oY (1.5) 

the subset of X of fully segregated pairs. Then, we study the ground states of <?gi ,32,00 when 
gi and (72 are both equal to zero, which results in an eigenvalue problem that leads to half- 
space geometry. 

Our first result is. 

Theorem A. Let (itg", Ugi)ngN be a sequence of minimizers of £gn over X, with = 
K, 52, 912) e [0,co? X M+ such that 

hm gi = gi, lim 52 = 92 and lim gi2 = oo. (1.6) 

n— >oo n— >oo n— >oo 

There exist a limiting pair {uoo,Voo) G Y such that, up to a subsequence, 

(i) {ugn,Vgn) converges to (^00,^00) weakly in Hy{M?) x Hy{R'^) and in C;°^(M^) x C;°^(M^). 

(ii) {uoo,Voo) minimizes <fgi,g2,oo overY and satisfies weakly the system 

—Au + Vu + gi\u\'^u = Xu in {|iioo| > 0} 
< -Av + Vv + g2\v\'^v = nv in {|voo| > 0} (1.7) 
u-v = in M^. 

Here A and /i are respectively the limits of the Lagrange multipliers A„ and associated 
with (ugn^Vgn) by (1.3). 

This is the analogue of the results of Wei and Weth in [33] for solutions of the system (1.3) 
in bounded domains and without trapping potential. We obtain the corresponding results by 
adapting the techniques in this papers to our setting. In [33] (see also the work of Chang, 
Lin, Lin and Lin in [11]), the authors prove the segregation in the strongly coupled case and 
the local uniform convergence to a limiting pair solving the system (1.7). In [14], Conti, Ter- 
racini and Verzini study the equivalent of the energy £g^,g2.oo defined over bounded domains 
and without trapping potential. They prove the existence of minimizers and their Lipschitz 
regularity, and give extremality conditions in the form of a system of subsolution of ellip- 
tic equations. We emphasize that we address the problem that actually corresponds to the 
physical situation. In the previous mentioned works, no mathematical results are presented 
about the symmetry breaking for solutions of system (1.3). 

The main result of this article, and the improvement with respect to the previous works, is 
to prove the symmetry breaking and to give an accurate description of the behavior of the 
ground states, in the strongly coupled and weakly interaction case: 

Theorem B. Let {ugn,Vgn) be a sequence of minimizers of £gn over X , with g^ = (51,(72 '512) 
such that 5f — )■ 0, ~^ '^'^^ 9i2 ~^ Then, up to a subsequence, Ugn and Vgn converge 
in C|^g(M^) respectively to e^^+w'^ and e^^^w~ , where 

^ a; I a: I 

w,^(x) = —= uj (x ■ v) e ~ , u £ S (1-8) 

is a second eigenfunction of the harmonic oscillator —A -|- a;^|xp in M^, and 0^, 0_ are real 
constants. 
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The limiting functions w'^ and w~ are not radial since they are supported in the half-spaces 
{x • 1/ > 0} and {x-v < 0}. Theorem B says that in compacts domains of M^, the supports of 
Ug and Vg approach these half-spaces when 5i — ?• 0, g'2 — ?■ and gi2 — ^ oo. The desired result 
about the symmetry breaking follows immediately: 

Corollary 1.1. There are positive constants go and G, such that if {ug,Vg) is a minimizer 
of £g in X with 

max{gi,g2} < 90 and gu > G , 
then Ug and Vg are not radially symmetric. 

Theorem B follows directly from the local uniform convergence stated in Theorem A, together 
with an accurate description of the segregated minimizers in the non interacting limit when 
91= 92 = 0, given by: 

Theorem C. Let (uo,vq) be a minimizer o/ ifo,o,oo over Y. Then, (i) uq is supported in 
H = {xi > 0} and vq in or rotated of them, and (ii) uq and vq are the positive and the 
negative parts of Wy, a second eigenfunction of the harmonic oscillator —A + in 
given by (1-8). 

We prove Theorem C using the results of Ehrhard in [15] about the extremality properties 
of half-spaces. Theorem C says that in the limit case when gi = g2 = 0, the nodal set of a 
segregated two-component BEC is a straight line. We emphasize that in the case when gi 
and g2 are small and positive, the exact configuration of both components is still unknown, 
but we expect the nodal set to be an infinite line with small curvature. We mention the work 
of Berestycki, Lin, Wei and Zhao in [3], where they study the profile of the components near 
the interphase using blow-up techniques. 

The study of the geometric nature of the support of the segregated minimizers of i?o,o,oo is 
an example of an optimal partition problem, or of a free boundary problem between several 
components. 

The optimal partition problem consists in finding a partition of C M*^ by k disjoint open 
sets uji, . . . ,ujk, that minimizes a function of the groundstate energies in each component. 
This is an open problem in general, and only few results are known about the exact con- 
figurations of these optimal partitions. In [9], Caffarelli and Lin consider the problem of 
minimizing Yli=i '^('^Oj where A(a;) is the first Dirichlet eigenvalue of the Laplacian on to. 
For 0, bounded, they show the existence of classical solutions as well as the regularity of the 
interfaces. In [20], Helffer, Hoffmann-Ostenhof and Terracini also consider the problem of 
minimizing maxj=i^...^fc /^(wj), where fJ,{uj) is the first Dirichlet eigenvalue of some Schrodinger 
operator on uj. They study the relation between nodal domains, spectral partitions and spec- 
tral properties in a bounded domain. An important question is whether a minimal partition 
is nodal, that is, consisting in the nodal domains of an eigenfunction of the Dirichlet realiza- 
tion of the Schrodinger operator in Q. This is the case for the second problem when k = 2 
and Q is bounded (see [14] and [19]). 

In Theorem C, we actually prove that the optimal partition, for the minimization of the 
addition of the energies of a system of two harmonic oscillators in M?, is composed of the 
supports of the positive and negative parts of its second eigenfunction. 
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The free boundary problem, which usually arises for competitive system when gi2 goes to 
infinity, has been studied by Cafi^arelli and Lin in [10], for the singularly perturbed elliptic 
system —Au + gi2u'^v = and —Av + gi2v'^u = 0. They prove the existence of a limiting pair, 
for which each component is a harmonic function in its domain of definition; and the C"'" 
regularity of the nodal line {u = v = 0}. They use Algrem's monotonicity formula, originally 
established in [2], which allows to prove Holder uniform estimates for the solutions. See [8], 
[13] and [27] for other works in this spirit. 

A further related problem is the study of (1.3) in the attractive case when gi and 52 are 
both negative. This is certainly a very different problem, since in this case the associated 
energy may not be bounded from below (see for example [23]) and no ground states exist. 
Symmetry results have been proved in this case, when there is neither trapping potential nor 
mass constraints. In [32], Wei and Weth prove that in a system defined over all M^, they are 
infinitely many non radial solutions, and in [34], that in a system defined in the unit ball, 
for every G N* there is a radially symmetric solution {u, v) such that u — v changes exactly 
k times of sign in the radial variable. We also mention the work of Liu in [24], where he 
does not prescribe the constraint of the L^-norm of each component being one, and therefore 
in the limit obtains that only one component remains. This is another form of segregation, 
since the model allows one component to disappear in the strong interacting case, when gi 
and g2 go to infinity. 

Sketch of proofs 

We now explain the outline of the paper and the main ideas of the proofs. 

In Section 2, we prove Theorem A and the segregation of the minimizing pairs {ug,Vg). We 
first remark that for all g £ and {u,v) G Y, £g{u,v) = <5g^^g2,oo('w, t"), so £g{ug,Vg) < 
£gi,g2,oo{u,v), and there are positive constants ci, C2 and C3 such that for every gi > and 
92 > 

£giug,Vg) <ci + C2gi + C3g2, (1.9) 

for all gi2 > 0. 

Sketch of the proof of Theorem A: Using (1.9), the energy of {ugn, Vgn) (and also the Lagrange 
multipliers Xg^ and /ig") is uniformly bounded. Hence, there exists (^00,^00) S X which is the 
weak limit of {ugn,Vgn) in Hy x Hy and the strong limit in !?■ x L?'. We prove that (uqo, Woo) 
minimizes £g^^g2,oo in Y using energy estimates. In order to prove the Cf^^ convergence, we 
follow the ideas Wei and Weth in [33], where they prove the equicontinuity of solutions of 
system (1.3) in bounded domains and without trapping potential. We show that the proof 
of Theorem 1.1(a) therein works in our setting. The reason of this, is that the proof consists 
in a rescaling of the solutions, which yields in a limit problem over all M^. Rescaling Ugn and 
Vg'n. identically, we get the same limit problem, so the equicontinuity holds, which gives the 
local uniform convergence. To show that (uoc'^oo) satisfies the system (1.7), we first prove 
some estimates for minimizing pairs {ug,Vg) of £g: there are C2 > and C3 > such that, 
jjiiglloo, ll^^glloo < C2 and [[Viiglloo, llVujjIloo < Csy/g^, for any g £ [0,co]^ x M+ with gu large 
enough. Using these, we show in Proposition 2.3 that there is C4 > 0, such that for every 
£ > and every > 1, \vg\ < Cig^2 ™ i^^^g'^^g > ^} l%l — ^^912 ™ {^^ff'^s > 
for gi2 large enough depending on e and rj. This implies that the functions 9i2\vgA'^Ug'^ and 
9i2\'Ugn\^Vgn converge weakly to zero respectively in {u^ > 0} and {voo > 0}, and the system 
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(1.7) is satisfied. 



Remark 1.2. We expect the derivatives of Uoo and v^o to be not continuous trougli the nodal 
hne {uoo = = 0}, and the system (1.7) to be solved only in the interior of the supports 
of Uoo and foo. In [22] and [25], we find simulations supporting this idea. In addition, in 
the case when gi = g2, we conjecture that the difference of the two components is a smooth 
function. Indeed, as we explain in Remark 3.6, in this case, we expect Uoo — t^oo to solve a 
homogeneous elliptic equation, which after standard elliptic regularity arguments, yields in 
smoothness for Uao — Voo ■ 

In Section 3, we study the properties of fully segregated BECs. We show that the minimizers 
of ^9i>92,oo in Y are locally Lipschitz continuous, and that the nodal set has empty interior. 
The local Lipschitz continuity is an important result because later in the proof of Theorem C, 
we need the limiting function to be locally Lipschitz continuous. The result about the nodal 
line of the segregated minimizers is used in the proof of Theorem C. We stress here that the 
space Y is not a manifold, so we cannot perform calculus of variations therein. We have to 
use other techniques to deal with minimizers of £g^^g2,oo in Y ^ in order to get a system of 
equations allowing us to study their local properties. We will explain these techniques, based 
on the works of Conti, Terracini and Verzini in [12] and [14]. 

Remark 1.3. All the results from Section 2 and Section 3 do not depend on the specific 
form of the trapping potential, and they remain true when considering a potential V with 
polynomial growth at infinity, i.e., a real function V such that 

V{x)>c\x\P iliM^\Br{Q) 
for some c > 0, r > and p > 0. 

In Section 4, we prove Theorem C and Theorem B. Theorem B follows directly from Theorem 
A and Theorem C(ii). The key ingredient in the proof of Theorem C is the extremal property 
of the half-spaces Ha,v = {x ■ v > a} with respect to the Gaussian-Rayleigh quotient 

Here a is a real number, v ^ and ^ is the Gaussian measure over M^. In [15], Ehrhard 
prove, roughly speaking, that if A(5') is the infimum of F over the functions vanishing outside 
of a Lebesgue measurable set 5, and Ha,u is a half-space with the same Gaussian measure as 
S, then 

K{S) > K{Ha,u) , 
and the equality holds only when S = Ha^v for some u ^ . 

Sketch of the proof of Theorem C: For the first assertion, we start noticing that the change 
of variables u{x) = \/27r a M(ax)e4l^l , whit a = (2a;)~^/^, gives 

^o,o,oo(n, ';;) = a; {F{u) + F{v) + 2) . 
Next, we first show in Lemma 4.2 that if {uq,vq) minimize fo,o,oo over Y, then 

F{uo) = HUo) and F(no) = A(Vo) 
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where Uq and Vq are respectively the supports of uq and vq. Using Ehr hard's result, we 
establish then that 

^0 = Ha,u and Vq = Hi,^yi 

for some a, 6 G and z/, i/' G 5^. Because ZYq ^-ncl Vq are disjoint, and because the nodal 
line has empty interior (see Proposition 3.4), we get that a = —b and u' = v. Hence 
Vo = H^j^. Then, we use a result of Beckner, Kenig and Pipher (see [7]) saying that the 
mapping a i— )■ h.{Ha^u) is convex, to prove that a is equal to zero. The definition of the change 
of variables gives then that suppuo = ifo,y and suppvQ = Hq^ for some G S^. Finally, 
using standard arguments for one-component BECs, we prove that system (1.7) is uniquely 
solved in half-spaces. The second assertion of Theorem C comes then easily, remarking that 
the positive and the negative parts of a second eigenfunction of the harmonic oscillator solve 
(1.7) in half-spaces. 

2 The segregation limit 

In this section we assume that gi and 52 are bounded with respect to 512, so there is cq > 
such that 

max{ ^1,52} < Co . (2.1) 

Moreover, with out loss of generality, we assume that Ug and Vg are real positive functions 
over all M^. We are allowed to do this after the following result, which is standard for one 
component BECs (see [1]): 

Lemma 2.1. Each component of a minimizing pair of 8g over X is, up to a complex multiplier 
of modulus one, a real positive smooth function. 

Proof. Every minimizing pair (ug,Vg) solves the system (1.3), so using standard elliptic 
regularity and the strong maximum principle, Ug and Vg are non vanishing smooth complex 
functions. Thus, there are smooth real functions ipi and ip2, such that Ug = \ug\e^^^ and 
Vg = \vg\e^'^^. The diamagnetic inequality imply that (|ng|, |ug|) is also a minimizer of £g over 
X. We have then the equality i5g(|ng|, |ug|) = £g{ug,Vg), which imply that ipi and (/?2 are 
constants, and hence the result. □ 

We start showing uniform estimates on £g and on Ug, Vg and its derivatives. We will use 
these to prove Theorem A and the segregation of Ug and Vg in Proposition 2.3. 

Estimates on minimizers 

Lemma 2.2. There are positive constants Cq, Ci, C2 and C3, such that if {ug,Vg) is a min- 
imizing pair of £g over X, with g = {gi, g2, 912) G [0, cq] x [0,co] x and Xg, fig are the 
associated Lagrange multipliers, then 



£g{Ug,Vg) < Cq 
0<Xg,Hg < Cl 

linllnoi 1 1 foil 00 ^ C2 



•^3 1 1 00 ) 
-^glloO) 



IVUolloo, llVUolloo < C3y/g^ 



(2.2) 
(2.3) 
(2.4) 
(2.5) 



for every gi2 > 1 . 
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Proof. Proof of (2.2): as we saw in (1.9), if {u,v) G Y then 



£g{Ug,Vg) < £g{u, v) = £gj,g2,oo{u, V) , 

so there are positive constant cq, ci and C2, not depending on gi, §2 or gi2, such that for 
every 51 > and 52 > 

£g{Ug,Vg) < CQ + ClQl + 02^2 = Cq , 

for aU 512 > 0. We get then (2.2) after (2.1). 

Proof of (2.3): Multiplying the first equation in (1.3) by Ug and then integrating over all 
we get that 



SO after the mass constraint Xg > 0, and after (2.2), Xg < A£g{ug,Vg) < ACq = Ci. The same 
argument is valid with Hg, which yields (2.3). 

Proof of (2.4): Consider x G and R > 0. Using (2.2), the mass constraint and the 
continuous embedding for p S [2, cxd), for every ball B = B2r{x) there is a positive 

constant C = C (p, R) such that 



\\ug\\LP(B) < C'\\ug\\Hi(B) < C'^2£g{ug,vg) + 1 < C'^^ICq + I. (2.6) 
After (1.3) we get 

— Aug < hg in B 

with hg = XgUg. Uslug (2.3) together with (2.6), there is C' = C" {q, R, gi, §2) > such that 

\\hg\\Li/'^(B) < C (2.7) 

for every q E [4, 00). 

Using a local estimate for subsolutions of elliptic equations (see Theorem 8.17 in [17]) 
there is C = C {R, p,q) > such that 

sup Ug < C'" {R^'Hu1,\\LnB)+R'~'Hhg\\L./HB)) ■ 
Br(x) ^ ^ 

Fixing i?, p and g, we derive from (2.6) and (2.7) that there is C2 > such that 

Ug{x) < C2 

for every gi2 > 0. The same argument is valid with Vg. Therefore, since Ug and Vg are 
positive, we get (2.4). 

Proof of (2.5): we first prove that Ug and Vg have polynomial decay at infinity. More precisely, 
we claim that for every a > there is > and > such that for all gi2 > 0, 

ugix) < ^ (2.8) 
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for all X eKa = M.^\ -B,.^(0). For a fixed, take = ^jiu^^iCx + ^JC'i^ ^o?u'^). A straight- 
forward calculation shows that 



\x\ 



fa{x) = C2 

is a supersolution of the first equation in (1.3), 

-AU + faiV + gif^ + gi2\v\^-Xg)>0 in 



while 



fo 



= C2>Ug 



Now define ip = fa — Ug and suppose that ^ is strictly negative somewhere in K^- Because 
fa and Ug are of class in Ka and go both to zero at infinity, ijj must have a local minimum 
xq in Ka ■ Ug{xo) > fa{xo) and D'^ijj^xo) is positive defined, so A'(/'(xo) > 0. Using this, and 
the fact that fa is a super solution of the first equation in (1.3) while Ug is a solution, we 
have 

{V{xo) + gifaixof + gi2Vg{xof - Xg)fa{xo) > {V{xo) + giUg{xof + gi2Vg{xof - Xg)ug{xo) . 

But our choice of Tq, implies in particular that V{xo) — > 0, so we get a contradiction 
with ^p{xo) < and the claim is proved. Remark that the previous claim, together with (2.4) 
imply that Vug is uniformly bounded in with respect to gi2- 

To finish the proof, let x £ M? and suppose that 512 > 1. For y G i?2(0) define u{y) = 



1/2 

u{x + 5i2 y) and v{x) = v{xo +5i2 u)- We have 

An(y) = gi2{y{x + gi2^y)u{y) + giu{yf - Xgu{y)} + u{y)tP'{y) , 

so after (2.1), (2.3), (2.4) and the claim there is a constant c > such that |Au(y)| < c for 
all y € -62(0) and gi2 > 1. Using a Holder estimate for the first derivative of u (see Theorem 
8.32 in [17]) there is a constant C > such that 

1 1 Vtt 1 1 LOO (5^(0)) < C{C2+C) . 

We get then the result with C3 = C(C2 + c) considering Vu (0). 

□ 

We now show the key ingredient in the Proof of Theorem A, the segregation of Ug and Vg. 
This is a generalization of Proposition 2.1 in [11] to positive solutions of (1.3) defined in all 
M? . We have also used in the proof some ideas from [33] . 

Proposition 2.3. Consider {ug,Vg) a sequence of minimizers of 8g over X and define 
Ue = {xeM^] miug[x)>e], = {x G ; mivg{x)>e]. 

3i2>0 3i2>0 

For any e > and rj > 1, there are Gq > and a positive constant C4 such that 

Vg < C4 g^2 o,nd Ug < C4 g^2 

for every gu > Gq. 
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Proof. Let i] > 1, e > and x G ZY^. Define the numbers 

ryi = fr? p e (0,^6-^2 



so Sg < tg for gi2 large enougli depending on rj and e, and note B = Bsg{x) 
Define also the function hg : (0, 00) — )• M by 



and notice that 



27rr h' Jr) = Ug ■ d^Ug ds . (2.9) 

JdBr{x) 

After (2.4) we have that < hg < C|. We claim the existence of ^g E {sg,tg) such that 
If not, we will get 



C| > Vt,)-^W>/;;=^* = ln(|^) 

^ / ln(r/i lngi2)-| lnffi2 \ 912^^^ 
1^ -pin J 

which contradicts the choice of p. 
We have then 

\Vug\'^ dx < / \Vug\'^ dx = / Ug-dijUgds— I Aug-Ugdx. 

IB JB^gix) JdB^gix) JB^gix) 

The first term of the right hand side can be estimated using (2.10), and the second one using 
(1.3), (2.3) and (2.4). We get 



l^\Vug\^dx<-^^+CiC^fg, 
so there is a positive C' such that 

r C' 

/ iVuJ'^dx < (2.11) 

Jb 111512 

for gi2 large enough. 

Using Theorem 7.17 in [17], we have that for p G (2,3) and 7 = 1 — |, there is Cp > such 
that 

OSCBUg < CpS2\\VUg\t-^^/P\\VUg\\l^r 
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Using (2.5), (2.11) and the definition of Sg, there is C" > such that 



_ J. _|„ A 

oscBUg < C rflg^2 ^ ''(111512)' 
= C" rfl{\ngi2f~p , 

so after (2.1), oscs — )• when gi2 — )• 00. This imphes that 



Ua> - in B 

for gi2 large enough. 

Using this last estimate, together with (1.3) (2.3) and (2.4), we get 

-Aug < -gi2^Vg in B 
in B 
in dB 



Vg>0 



(2.12) 



Vn<C2 



for gi2 large enough. Hence, Lemma 4.4 in [12] gives that exist a constant C > (not 
depending in e, or x), such that 



'912^ 



m 



Bf^{x) 
2 



Tacking C4 = C C2 , the definition of Sg gives 



Vg{x) < CAgi2 ) 

for gi2 large enough depending on 77 and e. The equivalent argument holds for Ug in V^, which 
yields Proposition 3.1. 

□ 

We have now all the tools to prove Theorem A. 

Proof of Theorem A. Let (ti„, Vn) = (ugn^Vgn) be a sequence of minimizing pairs of £gn in 
X with gf ^ gi, g2 ^ 92 and g'^^ ^ oo- 

(i) After (2.2) the sequences Un and f„ are bounded in Hy, so there exists (uoo'^oo) S 
X Hy with (up to a subsequence) 



Un 
Vn 



Uoo in Hy 
Voo in 



as n — 7- 00. The compact embedding Hy{M.'^) ^ L'^(M?) gives the strong convergence, so 
lliioolb = Iboolb = 1 and (Moo,foo) G 

On the one hand, we have that 

\\UnVn - '^^oo'^^oolll < H^inlbHWri - foolb + |b oo lb Ulin " ""oolb = o{l) , (2.13) 
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so (up to a subsequence) UnVn — )• ^^oo'^oo a-e- in M?. And on the other hand, after (2.4) and 
(2.8), UnVn < Cih G L^. The Lebesgue dominated convergence theorem together with (2.2) 
gives then 

||^toofoo||2 = hm llun^nlb < hm = . 
Hence, Uqo • t'oo = a.e. in and {uoo,Voo) S Y. 

In order to prove the C^^^ convergence, we follow directly the ideas of Wei and Weth in [33]. 
In Theorem 1.1(a) therein, they show that sequences of positive solutions of a class of com- 
petitive nonlinear elliptic systems in bounded domains are uniformly equicontinous. Their 
proof consists in a rescaling of the solutions and the domains, which yields in a limit problem 
over all . We will show that defining the rescaled functions of Ugn and Vgn , we get the same 
limit problem, so the equicontinuity holds. We recall that {un,Vn) satisfies the system (1.3) 
and that after Lemma 2.2, the sequence is uniformly bounded in Hy x Hy and in L°° x L°°. 

Following the proof of Theorem 1.1(a) in [33], if the sequence (ti„, f„) is not uniformly equicon- 
tinuous, there exists 6 > such that, without loss of generality, ii„ satisfies (up to a subse- 
quence) 

inf {[x — y| ; x, y G , |u„(x) — u„(y)| > 2(5} — as n — )• oo. 
Then, since the ti„ functions are positives, (2.8) implies that there are XmUn G such that 

Vn ■= \Xn - ?/n| as n oo, d„ := UniUn) > S and Un{Xn) >dn + S. 

Take ei = (1, 0) and choose An G 0(2) such that AnCi = r~^{yn — Xn)- We define the rescaled 
function Vi^n : M+ by 

Vl,n{x) = Un{Xn + TnAnV) and V2,nix) = Vn{Xn + TnAnV) . 

Then, vi^n solves in 



-Aui_„ = ll^nVl,n - rl f l,n 

Vl,n > 

Vl,n{^l) = dn > 6 

Vl,niO) > Wl,n(ei) + (5. 



(2.14) 



Here h,n{x) = {V{xn + r„^„,a;) + g^vj ^ix) - A„,), which after (2.8), (2.3) and (2.4) satis- 
fies 

h,nVi,n ^ in L°° as n — )• oo . (2.15) 

Moreover, after (2.2) fi^„ is uniformly bounded in H^, and Lemma 2.3 also applied for the 
sequence {vi^n,V2,n)- This last two properties, together (2.14), (2.15), implies that (fi,ni^^2,n) 
satisfies the same hypotheses as in the proof of Theorem 1.1(a) in [33]. Hence, we obtain 
the same limit problem when n goes to infinity, which following exactly the proof, yields a 
contradiction. The desired result then holds. 

(ii) For the first assertion, let (tt, v) be any pair in Y. Then 

£g'l,g^,oo{u,v) = £gn{u,v) > £gn{Un,Vn) > £g^ ,oo{Un, Vn) ■ (2.16) 
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Since the pair (u„, u„) satisfies the uniform bounds (2.8) and (2.4), the L?' convergence imphes 
the convergence. This, together with (2.16) and the weak lower semicontinuity of the Hy 
norm, gives 



n— >oo ^ ^ 



> \im.ml£g^gr^oo{Un,Vn) 



— '^91,92,00(^005^00); 

which imply the result since after (i), (1*00,^^00) £ Y- 

For the second assertion, let if he a, C°° function supported in CC {uoo > 0}. Multiplying 
the first equation on (1.3) by ip and then integrating, we get 

/ Vip ■ Vn„ + if {VUn + Qiui - XnUn) = - 5^2 ^^n ^^n ■ (2.17) 

JK JK 

Using the weak convergence of Un to Uoo, the left hand side of (2.17) tends to 

j Vif ■ VUoo + f (VUoo + 51 ■"L - '^^oo) 

with A the limit (up to a subsequence) of A„, which exists because of (2.3). 

After (i), n„, converges uniformly to Uqo in K. Hence, K CC {uoo > 2e} for some e > 0, and 
there exists N > such that 

K <Z \ inf > e} . 

n>N ^ 

Thus, Proposition 2.3, together with (2.4), implies that ^n^n converge uniformly to zero 
in so the right hand side of (2.17) tends to zero as n — ?■ 00. Hence, 

j Vv? • VUoo + f> {VUoo + giul^) = ^ j V^'oo • 

The same argument is valid with foo which yields the result. 

□ 



3 Properties of fully segregated two- component BECs 

In this section we prove some properties of fully segregated two-component BECs. The results 
in here will be used to prove the symmetry breaking in the limit case when 51 = 52 = 0. We 
start with the local Lipschitz continuity of minimizers of £g^^g2^oa in ^• 

Local Lipschitz continuity 

Proposition 3.1. If{u,v) is a nonnegative real minimizer Eg^^g^^oo overY, then u and v are 
locally Lipschitz continuous in . 

To prove this proposition, we first see in Lemma 3.3 that each component of a minimizing 
pair of £gi^g2,oo is a weak subsolution of a Gross-Pitaevskii equation, and that the difference of 
both components is a weak solution of a non homogeneous elliptic equation. We derive local 
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uniform estimates with respect to gi and g2 for the minimizing pairs. These estimates imply 
that the minimizers satisfy some extremahty conditions. We conclude using the following 
result of Conti, Terracini and Verzini ([14], Theorem 8): 

Theorem 3.2. (Conti- Terracini- Verzini) Let Q, be a bounded regular set of IS?, M > 
and wi,W2 G H^{^1) such that wi > 0, W2 > and wi ■ W2 = 0- If wi and W2 satisfy 

-Awi < M , -Aw2 < M and - M < -A{wi - W2) < M , 
then they are both Lipschitz continuous in the interior ofQ. 

Lemma 3.3. Let gi > and g2 > 0. If {u,v) is a nonnegative minimizer of Sg^^^g^^oo over Y 
then 

- Au + Vu + giu^ < Xu, - Av + Vv + g2V^ < fiv (3.1) 

and 

- A{u -v) + V{u -v) + giu^ - g2V^ = \u - jiv (3.2) 
in P'(M^), where A = ei(u) and fi = e2{v) with 



ei{w) = j \Vw\'^ + V\w\'^ + gi\w\ 



for i = 1,2. 

Proof. Proof of (3.1): Arguing by contradiction, suppose that 



J Vn • V0 + {Vu + giu^ - Xu) (j) > 



for some < (j) e C§° {M.'^) 



For t G (0, 1), define a new test function as : 

{u - t(py 



{wi,W2) 



{u - t(l))- 



Where n"*" = max(u, 0) and u = max(— ti,0). In the rest of the proof the o(-) notation will 
mean with respect to the t — t- limit. 

Since {{u - t(j))^ > 0} C {u > 0}, (n - t0)+ • t; = a.e. in M^, and {wi,W2) G Y. 
Using = (/"^)^ + (/~)^, ||u||2 = 1 and < (u — i(/))~ < tcf), we compute 



{u-t4>)m = i+y {[{n-m+?-u^) 

= I- j {2t(t)u + [{u-t(P)~f -t^<j)'^) 
= I- j 2t(j)u + o{t) , 



so 
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\\{u-H)ni 

The difference between the energies is then 



1 + / 2t(t)u + o{t) and — ^ ^,_^„. = 1 + / 4t(^u + o{t) , 

7 \\[u-t<j))+\\l J 



£g„g2,oo{wi,W2) - £g„g^,oo{u,v) = \ I (| V(n - - | Vnj^) + / t0n • / |V(n-t</.)+'2 



2 



+i / V ([(^ _ t^)+]2 -u') + I tc^u ■ I V{x)\iu - t<A)+P(3.3) 

1 

4 



We note that 



l(|V(n-i0)+p-|Vn|2) < |l(|V(n-t<^)|2-|Vnp) 



-t J V4>-Vu + o{t) , 



so (3.3) becomes 



£gi,g2,ooiwi,ui2) - £gT^,g2,ooiu,v) < -t j Vu-V(j) + Vu4> + giv?(j)- ei{{u-t(j))'^)u4> + o{t) . 

Using Lebesgue dominated convergence theorem, we have that ei{{u — t(j))~^) — )■ ei{u) when 
t — )■ 0, so 

■^91,92,00(^1, '(^2) - Sgug2,ooiu, v) < -t J Vu-V(l)+ {Vu + QiU^ - Xu) (j) + o{t) . 

Hence, for t small enough we get the contradiction 

£gi,g2,oo{wi,W2) - £g^^g^^oo{u,v) < 0. 

Using the same arguments with v and the inequalities in (3.1) are proved. 
Proof of (3.2): Define u = u — v and suppose that 

j Vu-V(j) + {Vu + giu^ - g2V^ - Xu + vv)<l) < 
for some 0<^eC^{M.^). 

For t G (0, 1), define a new test function as : 

, , / (n + t0)+ (u + tcf))' 

[Wi,W2) - 

As before, we compute 
1 



{u + t(P)+\\i 



{u + t(t))+\\2 ' \\{u + t(t>)-\ 
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1 + / 2t<j)v + o{t) and ,,. , = 1 + / 4t0t; + o{t) 



\\{u + t(t>yu 

Using that u • v = 0, the difference between the energies is 



+ \ I gi{\(.u + tcp)+\'-\u\') + g2{\{u + t 

-t {uei{{u + t(l))+) - V e2{{u + t(l))~ ))(/) + o{t) . 



Hence, 



£91,92,00(^^1,102) - Sgi,g2,ooiu,v) = t j Vu-V(j) + Vu(l) + (giu^ - g2V^)(l) 

- t (uei((u + t^)+)-ve2((u + t()))-))()) + o(t). 



Using the same argument as before, we see that ei((u + t(p)~^) — ei(u) = o(l) and e2((u + 
t(p)~) - e2(v) = 0(1), so 

£gug2,oo(wi,'W2) - £gj^,g2,oo(u,v) = t j V u ■ V <j) + (Vu + giu^ - g2V^) , (j) 

— tj (u X — V fi) (j) + o(t) . 

And again, for t small enough we get 

£gi,g2,oo(wi,W2) - £g^^g2,oo(u,v) < 0, 

a contradiction. 

We have proved the inequality 

V(n -v)-V(l) + V(u - v)(j) + (giu^ - g2V^)(j) - (An - i'oov)<l) > . 
Using the same arguments with v = v — u we get 

V(v - n) • V0 + V(v - u)4> + (g2V^ - giv?)(l) - (-\u + Voov)4> > 
for every < (/> G C°°(]R2). Equality (3.2) is then proved. 

□ 

Proof of Proposition 3.1: Let O be any bounded set of M^. After (3.1), u and v are 
respectively subsolutions of — An = An and —An = fiv. Arguing as in the proof of (2.4) 
n and v are uniformly bounded in with respect to gi and g2, and since they minimize 
£gi,g2,oo, A and /i are also uniformly bounded. There are then positive Mi,M2 = 0(gi,g2) 
such that 
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-Au < Ml and - Av < M2 in Q , 
so after (3.2) and the previous estimates, there are positive M^jM^ = 0{gi,g2) such that 

-A/3 < -A{u -v)< Mi ill ^ ■ 

Theorem 3.2 with Ad = max{Mi, M2, M3, iVf4} imphes then that u and v are Lipschitz con- 
tinuous in Q, and the result is proved. 

□ 

The nodal set 

We now prove that the nodal set has empty interior. We do this by following an idea of 
Chang et all. in [11]. The point is that if the nodal set has a ball B contained in it, then 
u can be extended to a solution u of an elliptic equation in suppti U B, but then, after the 
strong maximum principle and the mass constraint, u cannot vanish in the interior of its 
support. 

Proposition 3.4. If{u,v) is a nonnegative real minimizer o/'?gi,g2,oo overY, then the nodal 
set {x £M? ; u{x) = v(x) = 0} has no interior points. 

Proof. After Proposition 3.1, we know that u and v are locally Lipschitz functions, so 
Z// = {x G ; u{x) > 0} and V = {x G ; > 0} are open regular sets. Define 
U = \ V and suppose that the nodal line has an interior point. Then lA C. lA and the 
Lebesgue measure of lA is less than the Lebesgue measure of lA. Let w be the minimizer of 

over the functions w in Hq{U) such that Jqw"^ = 1. 

Is clear that {w,v) G Y, so after £g^^g^^ao{w,v) > £g^^g^^oo{u,v) we get £g^flfl{u) < £g-^^Qfi{w). 
This imply that u solves —Au + {V + gi\u\'^ — X)u = 0, with A defined as in Lemma 3.3. 
Therefore, by elliptic regularity, u is a function in U, and using the strong maximum 
principle, u = because it vanish in the interior of U. This contradicts the mass constraint, 
so the nodal line has no interior points. 

□ 

Remark 3.6. As we said in the introduction, we expect the derivatives of u and v to be not 
continuous trough the nodal line. In the case when gi = g2, we also expect the difference 
of the two components to be a smooth function. Indeed, after Proposition 3.3, if (u, v) is a 
nonnegative real minimizer of ^31,32, 00 over Y, then u — v solves the elliptic equation (3.2) in 
M?. When gi = g2, we expect the Lagrange multipliers A and /x to be equal, and equation 
(3.2) to be homogeneous. Standard elliptic regularity theory implies then that u — v is a C°° 
function over all M?. These properties are verified when gi and g2 are both equal to zero: in 
Theorem C(ii) we show that in this case, u and v are respectively the positive and negative 
parts of a second eigenfunction of the harmonic oscillator L inM?, so u — v G C°°(M?), X = n, 
and there is a jump of the derivatives of u and v trough the nodal line. 
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4 The non interacting limit 



In this section we study the minimizers {uo,vo) of i5o,o,oo over Y. We prove Theorem C, i.e., 
that the two components are supported in complementary half-spaces meeting at zero, and 
that they are respectively the positive and negative parts of a second eigenfunction of the 
harmonic oscillator — A + a;^|xp in M^. 

For a G M and v £ we define the half space 

Ha,u = {x eR^ ; X ■ u > a} 

and we write Ha = -ffa,{o,i) • The main idea in the proof of Theorem C is the extremal property 
of half-spaces with respect to the Gaussian-Rayleigh quotient 

Here fj, is the Gaussian measure in M^, which density with respect to the Lebesgue measure 
is given by 

dfj,{x) = —e~2^^^^dx. 

We remark that the invariance of the Gaussian measure with respect to rotations gives that 
for every v, v' G 



K{Hay) = A{Ha,u) = MH'-a,u) , (4.1) 

and that since M 9 a i— t- fi{Ha^u) is an increasing function, for every Lebesgue measurable set 
S there is a real a such that n{S) = fi{Ha^u)- 

For a non empty open S, we define J-{S) as the class of nonnegative non zero functions 
absolutely continuous on lines with support in S, and A by 



A(S) = inf F(f) . 

In [15] and [16], Ehrhard study isoperimetric inequalities in Gauss spaces and introduce the 
Gaussian symmetrization, a variant of classical symmetrizations used to solve isoperimetric 
problems, such as the principal frequency of a membrane or the torsional rigidity of a bar (see 
for example [6], [29] or [30]). In [16], Ehrhard prove that among all subsets with prescribed 
Gaussian measure, half-spaces have minimal A: 

Theorem 4.1. (Ehrhard) Let S be a non empty open subset of M? and a G M such that 
niS) = fj.{Ha). Then, 



A{S) > A{Ha) 

and the equality holds if and only if S = Ha^u for some u £ S^. Moreover, the infimum in 
A{Ha) is attained by some f G F{Ha)- 

We recall that the trapping potential is given by (1.2), so after (1.4) the energy is 
£o,o,oo{u,v) = ^ I ||Vn|2 + wVlVrt + i / \\Vv\^ +uj^\x\'^\v\^} , 
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and that the class of minimization Y is given by (1.5). 

To prove Theorem C, we first perform a change of variable in order to deal with the mini- 
mization problem in a diff'erent setting. For (u, v) define 

u{x) = \/27rati(ax) e^'"^'^ and v{x) = ^/2n av{ax) e^'^^'^^ (4.2) 
with a = (2a; )~^/^, and 

Y = {{u,v)-{u,v)eY]. 

A direct computation gives 

^o,o,oo(^i, v)=uj {F{u) + F{v) + 2) . (4.3) 

Lemma 4.2. A nonnegative pair {u,v) minimizes Sq^q^oq overY if and only if{u,v) minimize 
F{u) + F{v) over Y . Moreover, in this case u minimizes F over J~{U) and v minimizes F 
over J-{y ), where lA = suppu and V = suppv. 

Proof. The first assertion is immediate after (4.3). If {u,v) is a nonnegative minimizer of 
^0,0, cx> over Y, then after Proposition 3.1 u is in C|^^(M^), so it is absolutely continuous on 
lines, and because of the mass constraint u is not identically zero. After (4.2) the same prop- 
erties hold for u, so u G F{U). For every w G -^(^) with finite Gaussian- Rayleigh quotient, 
{^/\\w\\2,v) G Y. Thus, (4.3) gives F{u) < F{w). The same argument is valid with v, so the 
second assertion is proved. 

We are now able to prove Theorem C. For the first part, we use the ideas of Beckner, Kenig 
and Pipher in Section 2.4 of [7]. 

Proof of Theorem C(i). To lighten the notation we write u = uq, U = suppu, lA = suppu 
and the analogous for vq. The diamagnetic inequality imply that (|ii|, |i;|) is also a minimizer 
of ^^0,0,00 over Y , so we suppose, without loss of generality, that u and v are nonnegative real 
functions. 

Step 1: U = Ha^u o-nd V = H^y ■ Suppose that or V is not a half-space. Then, after Lemma 
4.2 and Theorem 4.1 there are real numbers a, b such that /i(^) = fi{Ha), //(V) = fJ'iHb) and 

F{u) + F{v) = K{U) + A(V) > K{Ha) + A{Hb) . (4.4) 
We claim that a -|- 6 > 0. First, since n{U n V) = we have 

li{Ha) + fi{Hk) < 1 , (4.5) 

which imply that a and b cannot be both negative. We suppose then without loss of generality 
that 

b<0<a and a + b <0. (4.6) 

The fact that fJ.{Ha) + ^{Hta) = 1, together with (4.5), implies that iJ-{Hb) < fi{H^a), which 
contradicts (4.6). The claim is then proved. 

The inequality a -|- 6 > imply that Ha n H'^^^ = 0. Hence, for every pair {wi,W2) G 
J='{Ha) X F{Hif^), ("'i/ll«'ill2,'^2/ll«'2|j2) G Y. After (4.3) and (4.4) we obtain 
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F{wi) + F{w2) > F{u) + F{v) > A{Ha) + A{Hb) 



Minimizing F in the previous inequality with respect to wi G F{Ha) and with respect to W2 G 
J-{H'if^), and considering (4.1), we obtain the contradiction A{Ha) + A{Hi,) > A{Ha) + A{Hi)), 
so Step 1 is proved. 

Step 2: a = —b and u = v' . Since ■u • "D = 0, fM{Ha^u H H^^y) = 0, which imply that the 
boundaries of Ha^u and H^{,y must be parallel, i.e., u = v' and a > —b. Moreover, if o > 6, 
then the nodal set {u = v = 0} has a non empty interior, which after (4.2) contradicts 
Proposition 3.4. We have shown that 



U = Ha,u and V = H, 
for some a S ffi, S^. 

Step 3: a = 0. First, the monotonicity of the first eigenvalue of the Dirichlet problem with 
respect to the domain gives that a i— t- A[Ha) is an increasing function. Moreover, the same 
argument of Theorem 2.4.5 in [7] (see also Theorem 6.2 in [5]) gives that it is a convex 
function. Considering (4.1) we derive 

A{Ha) + A{Hl) ^ A{Ha)+A{H^a) ^ ^^^^^ ^ A(ffo) + A(gg) _ 

Suppose now that a 7^ and consider {wi,W2) G F{Hq) x F{Hq). The same argument used 
in Step 1, together with (4.7), gives 

F{wi) + F{w2) > A{Ha) + A(i^S) > A{Ho) + A{HS) . 

After Theorem 4.1, the infima in A{Hq) and A{Hq) are attained, so again, minimizing F in 
the previous inequality with respect to wi G F{Hq) and with respect to W2 G F{Hq), we 
obtain the contradiction A{Hq) + A{Hq) > A{Ho) + A{Hq). We have proved that 



U = Ho^^ and V = 



0,u ' 



which considering (4.2), gives 



U = Ho^u and V = ifg,^ 

for some v £ S^. 



□ 



We now give the proof of the second part of Theorem C, this is, that every minimizing pair 
of "^^0,0,00 over Y is of the form -v/2 (e*^+ w^, e*^" w~) for some 1/ G 5^ and 0^,6^ G M. 

The second eigenvalue of the harmonic oscillator —A + uj'^\x\'^ over has multiplicity 2. An 
orthogonal base, with respect to the product, of the associated spectral space is given by 

rii{x) = xi e 2 and mi^) = Cuj X2 e 2 

with Cuj = Y^V^"^- Every function in the spectral space, with norm equal to 1, is the of 
the form 
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Wu{x) = {x ■ ly) e 2 (4.8) 

for some z/ G 5^. 

Proof of Theorem C(ii). We write uq = u, = v and Hqu = H^. After Theorem A(ii) 
and Theorem C(i), there \s u ^ S"^ such that u solves 

— An + cj2|xpn = \u in Hy 

u = in dHy (4.9) 

ky = 1 

with ^ = fjj |Vup + uj'^\x\'^u'^ , and v solves the equivalent problem in H^. 

The same argument as in Lemma 2.1 imply that u = e*^+|n| in H,y and v = e^^~\v\ in H,y, 
with \u\ and \v\ positive and 6+,6- G M. We claim now that there is uniqueness for the 
modulus of u in problem (4.9). Suppose that there are two positive solutions ui and U2 of 
(4.9) respectively with Ai and A2. Suppose that Ai < A2 and define h = ui/u2 in H^. We 
have that 

V{ulVh) = {X2- Xi)ulh, (4.10) 
and using the mass constraint that 



ulh{h-l) = ]^l ul{h-lf 



Multiplying (4.10) by /i — 1, then integrating over and performing an integration by parts 
we find 



ul\Vh\'' + \{\2-\l)ul{h-lf = 



which imply that h = cte. The mass constraints imply then that h = 1, so the claim is 
proved. The equivalent result is valid for v in i/^. 

Finally, a direct computation shows that 

Uy{x) = V 2 Cuj {x ■ v) e ~ and Vi,{x) = v2cu) {x ■ u)" e ~ 

are positive solutions of problem (4.9) respectively in Hy and H^, sou = e^^+u^ and v = e^^~Vy 
for some u £ and 9+,6- £ M. 

□ 

We finish this article by given the proof of Theorem B. 

Proof of Theorem B. Let {ugn, Vgn) be a sequence of minimizing pairs of £gn, with (7" — )■ 0, 
^2 ~^ and — ^ CO- After Theorem A, Ugn (respectively Vgn) converges locally uniformly 
to uq (respectively uq), where (uc^o) minimizes £^0,0,00 over Y. Theorem C(ii) gives then 
that uq = e^^+w'^ and vq = e^^-w~ for some u e and 6+,9- G M. 

□ 

The author would like to thank C. Kenig, J. Wei, A. Aftalion and S. Terracini for their help 
and suggestions, and P. Mason for sharing with me his physical insights. 
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